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ABSTRACT 

Context. Evolutionary models of massive stars predict a surface enrichment of Nitrogen, due to rotational mixing. Recent studies 
within the VLT-FLAMES survey of massive stars have challenged (part of) these predictions. Such systematic determinations of 
Nitrogen abundances, however, have been mostly performed only for cooler (B-type) objects. 

Aims. This is the first paper in a series dealing with optical Nitrogen spectroscopy of O-type stars, aiming at the analysis of Nitrogen 
abundances for stellar samples of significant size, to place further constraints on the early evolution of massive stars. Here we con- 
centrate on the formation of the optical N m lines at /M4634 - 4640 - 4642 that are fundamental for the definition of the different 
morphological T -classes. 

Methods. We implement a new Nitrogen model atom into the NLTE atmosphere/spectrum synthesis code fastwind, and compare 
the resulting optical Nm spectra with other predictions, mostly from the seminal work by Mihalas & Hummer (1973, ApJ 179, 827, 
'MH'), and from the alternative code cmfgen. 

Results. Using similar model atmospheres as MH (not blanketed and wind-free), we are able to reproduce their results, in particular 
the optical triplet emission lines. According to MH, these should be strongly related to dielectronic recombination and the drain by 
certain two-electron transitions. However, using realistic, fully line-blanketed atmospheres at solar abundances, the key role of the 
dielectronic recombinations controlling these emission features is superseded - for O-star conditions - by the strength of the stellar 
wind and metallicity. Thus, in the case of wind-free (weak wind) models, the resulting lower ionizing EUV-fluxes severely suppress 
the emission. As the mass loss rate is increased, pumping through the N in resonance line(s) in the presence of a near-photospheric 
velocity field (i.e., the Swings-mechanism) results in a net optical triplet line emission. A comparison with results from cmfgen is 
mostly satisfactory, except for the range 30,000 K < r cff < 35,000 K, where cmfgen triggers the triplet emission at lower T sB than 
fastwind. This effect could be traced down to line overlap effects between the Nm and Om resonance lines that cannot be simulated 
by fastwind so far, due to the lack of a detailed O in model atom. 

Conclusions. Since the efficiency of dielectronic recombination and 'two electron drain' strongly depends on the degree of line- 
blanketing/-blocking, we predict the emission to become stronger in a metal-poor environment, though lower wind-strengths and 
Nitrogen abundances might counteract this effect. Weak winded stars (if existent in the decisive parameter range) should display less 
triplet emission than their counterparts with 'normal' winds. 

Key words, stars: early-type - line: formation - stars: atmospheres - stars: winds, outflows 

1 . Introduction studies provide the first statistically significant abundance mea- 

surements of Galactic, LMC, and SMC B-type stars, covering a 
One of the key aspects of massive star evolution is rota- broad range of rotational velocities, 
tional mixing and its impact. Evolutionary models including 

rotation (e.g.. iHeger & LangeJ l200Cl iMevnet & Maedel [2000; For the Galactic case, the mean B-star Nitrogen abundance 



?S ; IBrottetal . 201 la)~predict a surface enrichment of Nitrogen with as derived by | Hunter et al. | <|2009ab is in quite good agreement 

b ■ an associated Carbon depletion during the main sequence evo- wlth the corresponding baseline abundance. For the Magellanic 

■ ■ ' lutionQ The faster a star rotates, the more mixing will occur, Clouds stars > however, the derived Nitrogen abundances show 

and the larger the Nitrogen surface abundance that should be ob- clea rly the presence of an enrichment, where this enrichment is 

served more extreme in the SMC than in the LMCo Theoretical consid- 

" Several studies dHunteretalJ 120081 l2009fl IBrottetal] erations have ma J or difficultie s in explaining several aspects of 

l20TTbh have challenged the predicted effects of rotational mix- the accumulated results: Within the population of (LMC) core- 

ing on the basis of observations performed within the VLT- Hydrogen burning objects, both unennched fast rotators and 

FLAMES survey of massive stars (lEvans et all I2006T) . These m 8 hl y enrlched slow rotators have been found ' ln contradiction 

to standard theory, as well as slowly rotating, highly enriched B- 

* Appendix A, B, and C are only available in electronic form at supergiants (see below). Taken together, these results imply that 

http-//www edpsciences org standard rotational mixing might be not as dominant as usually 

1 due to a rapid achievement of the CN equilibrium, whilst the 
Oxygen depletion implied by the full CNO equilibrium is onl y found in 

rapidly rotating and more massive stars at later stages, e.g. jBrott et alj 2 Baseline abundances for all three environ ments from Hn regions 

1201 lal and unevolved B-stars, see Hu nter et alj d2007l) . 



1 



J.G. Rivero Gonzalez, J. Puis, & F. Najarro: N in emission line formation revisited 



quoted, and/or that other e nrichment processes might be present 
as well (lBrottetalJl2011bh . 

Interestingly, there exist only few rapidly rotating B- 
supergiants, since t here is a steep drop o f rotation r ates below 
T eS « 20 kK (e.g.. lHowarfhetal.lll997l) . Recently, IVink et all 
(2010) tried to explain this finding based on two alterna- 
tive scenarios. In the first scenario, the low rotation rates 
of B-supergiants are suggested to be caused by braking due 
to an increased mass loss for T e tf < 25 kK, related to the 
so-ca lled bi-stability jump (Pauldrach & Pulslll990t IVink et ail 
2000). Si nce the reality of such an increased mas s loss is still 
debated dMarkova & Pulsl 120081: iPuls et alJ l2010h . IVink et al l 
discuss an alternative scenario where the slowly rotat- 
ing B-supergiants might form an entirely separate, non core 
hydrogen-burning population. E.g., they might be products of 
binary evolution (though this is not generally expected to lead to 
slowly rotating stars), or they might be post-RSG or blue-loop 
stars. 

Support of this second scenario is the finding that the ma- 
jority of the cooler (LMC) objects is strongly Nitrogen-enriched 
(see above), and Vink et al. argue that "although rotating mod- 
els can in principle account for large N abundances, the fact 
that such a large number of the cooler objects is found to be 
N enriched suggests an evolved nature for these stars." A care- 
ful Nitrogen analysis of their (early) progenitors, the O-type 
stars, will certainly help to further constrain these ideas and 
present massive star evolution in general. Note that one of the 
scientific drivers o f the current VLT-FLAMES Tarantula survey 
(lEvans et al.ll201ll) is just such an analysis of an unprecedented 
sample of 'normal' O-stars and emission-line stars. 

Until to date, however, Nitrogen abundances have been sys- 
tematically derived only for the cooler subset of the previ- 
ous VLT-FLAMES survey, by means of analyzing Nn alone, 
whereas corresponding results are missing for the most mas- 
sive and hottest stars. More generally, when inspecting the avail- 
able literature for massive stars, one realizes that metallic abun- 
dances, in particular of the key element Nitrogen, h ave been de- 
rived only for a small number of O-type stars (e.g., Boure tetaT] 
l2003HHillier et alJl2003HWalborn et alJl2004l:lHeap et al.l l2006). 
The simple reason is that they are difficult to determine, since 
the formation of N ra/N iv lines (and lines from similar ions of 
C and O) is problematic due to the impact of various processes 
that are absent or negligible at cooler spectral types. 

One might argue that the determination of Nitrogen and 
other metallic abundances of hotter stars could or even should 
be performed via UV wind-lines, since these are clearly visi- 
ble as long as the wind-strength is not too low, and the line- 
formation is less complex and less dependent on accurate atomic 
models than for photospheric lines connecting intermediate or 
even high-lying levels. Note, however, that the results of such 
analyses strongly depend on the assumptions regardi ng and the 
treatm ent of wind X-ray emission and wind clumping ( Pui s et al] 
2008 and references therein). A careful photospheric analysis, 
on the other hand, remains rather unaffected by such problems 
as long as X-ray emission and clumping do not start (very) close 
to the photosphere, and we will follow the latter approach, con- 
centrating on optical lines. 

This paper is the first in a series of upcoming publications 
dealing with Nitrogen spectroscopy of O-type stars. The major 
objective of this project is the analysis of optical spectra from 
stellar samples of significant size in different environments, to 
derive the corresponding Nitrogen abundances which are key to 
our understanding of the early evolution of massive stars. In the 
present study, we will concentrate on the formation of the op- 



tical N in emission lines at /t/14634 - 4640 - 4642, which are 
fundamental for the definition of the different morphological T- 
classes. During our implementation of Nitro gen into th e NLTE- 
atmosphere/line formation code fastwind (Pui s et afll2005l) it 
turned out that the c anonical explanation in ter ms of dielec- 
tronic recombination (Mihalas & Hummer! 1 19731) no longer or 
only partly applies when modern atmosphere codes including 
line-blocking/blanketing and winds are used to synthesize the 
N m spectrum. Since the f-features are observed in the majority 
of O-stars and strongly dependent on the Nitrogen abundance, a 
thorough re-investigation of their formation process is required, 
in order to avoid wrong conclusions. 

This paper is organized as follows: At first (Sect. [2]) 
we recapitulate previous explanations for the Nm triplet 
emission, in particular the standard picture as provided by 
Miha las & Hummer] d!973l) . In Sect. [3] we discuss the dielec- 
tronic recombination process and its implementation into fast- 
wind. Sect. 3]presents our new N m model ion, and in Sect. [5] we 
investigate the dependence of the triplet emission on various pro- 
cesses. We compare our resu lts with corresponding ones from 
the alternative code cmfgen (iHillier & Miller! Il998h in Sect [6] 
and discuss the impact of coupling with Om via correspond- 
ing resonance lines in Sect. [7] The dependence of the emis- 
sion strength on specific parameters is discussed in Sect. [8] and 
Sect. |9]provides our summary and conclusions. 

In the next paper of this series ('Paper IF), we will present 
our complete Nitrogen model atom, and perform a Nitrogen 
abundance analysis for the LMC O-stars from the previous VLT- 
FLAMES survey of massive stars. 



2. Nm emission lines from O-stars - status quo 

The presence of emission in the N m triplet at AA4-634 - 4640 - 
4642 A, in combination with the behaviour of Hen /14686, 
is used for cl assification purposes ('f -features, see Walborn 
fl97Tl ISota et"aIll20Tlh . and to discriminate O-stars with such 
lin e emission from pu r e ab sorption-line objects. As stressed 
bv lBruccato & Mihalas ( 1971L hereafter BM71), these emission 
lines originate in t he stellar photosp here and not in an "exte- 
rior shell" (see also iHeap et alJl2006l for more recent work). In 
this case, the most plausible explanation is by invoking NLTE 
effects. In NLTE, line emission occurs when the corresponding 
source function at formation depths is larger than the continuum 
intensity at transition frequency vo. Such a large source function 
becomes possible if the upper level of the transition is (consider- 
ably) overpopulated with respect to the lower one, i.e., if b u > b\ 
(with b u and b\ the NLTE departure coefficients of the upper and 
lower level, respectively). Note that both coefficients can lie be- 
low unity. 

For the N m triplet produced by the 3p - 3d transitions (see 
Fig.Q]), such a mechanism should result in a (relative) overpopu- 
lation of the upper level, 3d. In the early work about Of stars, 
the fluorescence mechanism developed by iBowenl (1 19351) has 
been suggested to trigger such an overpopulation . Many authors 
([Swings & Sti-uve 1940; Swingsll948tlOkdl954 argued against 
the relevance of the Bowen mechanism in Of stars, because of 
the lack of Om emission lines at /L13340, 3444, 3759 A which 
are connected to the involved levels. 

An alternative suggestion is due to lSwingsl (Q948) and relies 
on an intense continuum that may directly pump the resonance 
transition 2p — » 3d (thus producing the required overpopulation, 
^3d > ^3p) as well as 2p — > 3s, while the transition 2p - 3p is 
radiatively forbidden. The implied overpopulation of level 3s 
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Fig. 1. Grotrian diagram displaying the transitions involved in 
the N in emission lines problem. The horizontal line marks the 
N in ionization threshold. The N in /L14634 - 4640 - 4642 triplet 
is formed by the transitions 3d — > 3p, while the absorption lines 
at /1/14097-4103 are due to the transitions 3p — > 3s. An effec- 
tive drain of 3p is provided by the 'two electron transitions' 3p 
— > 2p 2 ( 2 S, 2 P, 2 D). The levels above the ionization limit are the 
autoionizing levels that feed level 3d via dielectronic recombina- 
tion (Sect. [XT). Cascade processes (4f — » 3d) can overpopulate 
the 3d state as well. The Swings mechanism involves the reso- 
nance transitions 2p — > 3s and 2p — > 3d. Note that the (energetic) 
positions of the autoionizing niveaus have been shifted upwards 
for clarity. 



due to pumping may then explain why the transitions 3s - 3p 
at /L14097-4103 are always in absorption {bs s > £>3 P ). 

Without such pumping of the 3 s level, an auxiliary draining 
mechanism for the 3p level is needed, since otherwise an over- 
population, Z?3 P > £>3 S , might occur due to cascade processes, 
implying the presence of emission at /1/14097-4103, which is 
not observed. Indeed, the anomalous 'two electron transitions' 
2p 2 ( 2 S, 2 P, 2 D) - 3p with transition probabilities comparable 
to the 3p — > 3 s one electron transition have been identified by 
iNikitin & Yakubovskiil (Tl963) as potentially important draining 
processes. Calculations by BM71 show that the presence of these 
draining processes is sufficient both to ensure the overpopulation 
of 3d (relative to 3p) and to prevent the overpopulation of 3p rel- 
ative to 3s. Thus, these 'two electron transitions' play a key role. 

BM7 1 also noted a problem for the Swings mechanism when 
applied to realistic conditions. The 2p — > 3s and 2p — > 3d reso- 
nance lines are expected to be much more opaque than 3p — » 3 s 
and 3d — > 3p, and should be consequently in detailed balance in 
the line forming region. That would mean no pumping and no 
overpopulation of 3d relative to 3p (but see Sect. 15.31 ). The prob- 
lem became (preliminary) solved when BM7 1 suggested a third 
potential mechanism. They realized the existence of a large num- 
ber of autoionizing levels that either connect directly to the 3d 
state or to levels that can cascade downwards to 3d. Hence, the 
latter state may become strongly overpopulated by dielectronic 
recombination ('DR', see Sect.0. 

The most influential analysis of the Nm emission lines 
proble m until now has been carried out by Mi halas & Hummer] 
(I1973L 'MH'), building on the work by BM71. They used static, 
plane-parallel models trying to explain the effect for 0((f)) and 
O(f) stars. As a final result, they were able to reproduce the Nm 
triplet emission at the observed temperatures and gravities in 



parallel with absorption at AA4097 - 4103, by overpopulating 
level 3d primarily via dielectronic recombination. The subse- 
quent 3d — > 3p cascade produces the emission. The strong drain 
3p — > 2p 2 via 'two electron transitions' enhances the overpopu- 
lation of 3d relative to 3p by depopulating 3p and prevents emis- 
sion in the 3p — > 3s lines. Until to-date, dielectronic recombi- 
nation is the canonical explanation for the formation of the f- 
features. 



3. Dielectronic recombination 

If two electrons are excited within a complex atom/ion with sev- 
eral electrons, they can give rise to states with energies both be- 
low and above the ionization potential. States above the ioniza- 
tion limit, under certain selection rules, may preferentially au- 
toionize to the ground state of the ion plus a free electron. 

Thus the ionization from an initial bound state A(i) to an 
ionized final state A + (f) can occur either directly, or by a (dif- 
ferent) transition from the initial bound state A(i) to an interme- 
diate doubly excited state A above the ionization potential that 
finally autoionizes to A + (f) . In this respect, Pho to-Excitation of 
the Core (PEC) resonances dYan & Seatonl 19871) are particularly 
strong, because they correspond to a single electron transition in 
which the outer electron is a spectator and does not change. 

The inverse process is also possible, if an ion collides with 
an electron of sufficient energy, leading to a doubly excited 
state. Generally, the compound state will immediately autoionize 
again (large autoionization probabilities, A a ~ 10 13 - 10 14 s _1 ). 
In some cases, however, a stabilizing transition occurs, in which 
one of the excited electrons, usually the one in the lower state, 
radiatively decays to the lowest available quantum state. This 
process is the dielectronic recombination, and can be summa- 
rized as the capture of an electron by the target leading to an 
intermediate doubly excited state that stabilises by emitting a 
photon rather than an electron. 

3.1. The Nm emission triplet at AA4634 - 4640 - 4642 

As pointed out by Mihaiai dl97ll) . there are two important au- 
toionizing series in the Nm ion, of the form 2s2p( 1 P°)«Z and 
2s2p( 3 P°)n/, along with few bound double excitation states with 
similar configuration. Since states of the form 2s2pn/ are directly 
coupled to 2s 2 «/ states, they are of great importance. In the fol- 
lowing, we consider only the states from the singlet series, be- 
cause the transitions from the triplet series to 2s 2 «/ states are not 
electric dipole allowed transitions in LS coupling. 

The singlet series comprises only two bound configurations, 
2s2p3s and 2s2p3p, whilst the 2s2p3d configuration lies only 
1 .6 eV above the ionization potential and is of major importance, 
since the low position in the continuum produces strong dielec- 
tronic recombination, 

2s 2 ('S) (N iv) + e" -» 2s2p( 1 P°)3d( 2 P° , 2 F°). (1) 

Note that the core electron transition (see above), 
2s 2 ('S) — > 2s2p('P ), is equivalent to the resonance transi- 
tion in Niv. Thus, any 2s 2 ('S)nZ — > 2s2p( 1 P°)«/ transition 
gives rise to strong, broad resonances (see Fig. |2), unless the 
2s2p( 1 P°)«/ state is truly bound. 

The autoionizing states (Eq. [TJ can stabilize via two alterna- 
tive routes, either 

2s2p('P°)3d( 2 P , 2 F°) -> 2s2p( 1 P°)3p( 2 S , 2 P , 2 D) + hv, (2) 
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Fig. 2. Comparison of the 'raw' cross-section from the 
OPACITY project (black), and the smoothed one (grey/blue), 
for the 3d 2 D state of N m. Note the numerous complex reso- 
nances. An example for a PEC resonance, being much broader 
than the usual Rydberg resonances, is marked by the vertical 
arr ow. Dashed: c orresponding resonance-free data in terms of 
the lSeatonl d!958f) approximation from WM-basic. 



grid required by the resonances in the 'raw' data, which would 
increase the computational effort considerably. To circumvent 
the probl em we implemented a method that is also used in 
WM-b asic (fPauldrach et al.l200Th an d in cmfgen (Hillier & Miller 
119981) (see also lBautista et al.lll998h . The raw OPACITY project 
cross-sections are smoothed (Fig. [2]l to adapt them to the stan- 
dard continuum frequency grid used within the code, which 
has a typical resolution of a couple of hundred kms" 1 . To this 
end, the data are convolved with a Gaussian profile of typically 
3,000 km s _1 width, via a Fast Fourier Transform. This convolu- 
tion ensures that the area under the original data remains con- 
served, and that all resonances are treated with sufficient accu- 
racy. By means of this approach, the ionization and recombina- 
tion rates should be accurately represented, except for recombi- 
nation rates at very low temperatures, where the recombination 
coeffi cient is quite sensitive to the exact location of the resonance 
(see Hillier & Millejl998|). 

3.3. Implicit vs. explicit method 

Though under similar assumptions both methods achieve similar 
results, there are certain advantages a nd disadvantages that are 
summarized in the following (see also lHillier & Milled! 19981) . 



that end in the doubly excited bound configuration mentioned 
above, or the one that might overpopulate the upper level of the 
transition 3d — > 3p, 

2s2p( I P°)3d( 2 P° , 2 F°) -> 2s 2 3d( 2 D) + hv, (3) 

thus playing a (potential) key role in the N m triplet emission. 
The latter route is displayed in the upper part of Fig. Q] 

3.2. Implementation into fastwind 

Dielectronic recombination and its inverse process have been 
implemented into fastwind in two different ways, to allow us 
to use different data sets. Specifically, we implemented 

(i) an implicit method where the contribution of the dielec- 
tronic recombination is already included in the photoioniza- 
tion cross-sections. This method needs to be used for data 
that has been calculated, e.g., within the OPACITY project 
dCunto & Mendozall 19921) and 

(ii) an explicit method using resonance-free photoionization 
cross-sections in combination with explicitly included stabi- 
lizing transitions from autoionizing levels. This method will 
be used when we have information regarding transition fre- 
quencies and strengths of the stabilizing transitions. 

Later on we discuss the advantages and disadvantages of both 
methods, whereas in Appendix lA.ll we provide some details on 
the explicit method, and in Appendix IA. 21 we show that the im- 
plicit and explicit method are consistent as long as the autoioniz- 
ing states are in LTE, as often the case (or frequently assumed). 

Photoionization cross-sections from the OPACITY project 
data include the contributions of dielectronic recombination and 
will be used within the implicit method. All these cross-sections 
display complex resonances (where the largest and widest ones 
are the PEC resonances), which somewhat complicate the imple- 
mentation of this method. Since some of the resonances are quite 
narrow, care must be taken when sampling the cross-sections. 
If performing a straightforward calculation, the radiative trans- 
fer would need to be solved at each point of the fine-frequency 



The implicit method has the advantage that for states that can 
autoionize in LS coupling their contribution to dielectronic re- 
combination (and the inverse process) is already included within 
the photoionization data. There is no need to look for both the 
important autoionizing series to each level and the oscillator 
strengths of the stabilizing transitions. On the other hand, there 
are also disadvantages. Narrow resonances are not always well 
resolved, and, if the resonance is strong, the dielectronic recom- 
bination rate from such a resonance could become erroneous. 
Besides, the positions of the resonances are only approximate. 
If line coincidences are imp ortant, this can affect t he transfer 
and the corresponding rates. iHUlier & Milled ([T998) suggested 
to avoid the implicit approach for those transitions where dielec- 
tronic recombination is an important mechanism. In our new N m 
model ion (Sect. [4]i we have followed their advice. Finally, di- 
electronic recombination rates calculated by the implicit method 
are inevitably based on the assumption that the autoionizing lev- 
els are in LTE with respect to the ground state of the next higher 
ion0 In the rare case that this were no longer true, the implicit 
method cannot be used. Then, the autoionizing levels need to 
be included into the model atom, and all transitions need to be 
treated explicitly. 

The explicit method has the theoretical advantage that res- 
onances can be inserted at the correct wavelength if known. 
(Though one has to admit that resonance positions will never 
be very accurate. For the PECs this does not matter though.) 
Often, however, the profile functions for the resonances are dif- 
ficult to obtain (which, on the other hand, are included in the 
implicit approach). The width of these Fano profiles is set by the 
autoionization coefficients, which, frequently, are not available. 
To overcome this problem, we follow the approach by MH and 
assume the resonance to be wide (which is true for the most im- 
portant PEC-resonances), such that we can use the mean inten- 
sity instead of the scattering integral when calculating the rates, 



3 This will almost always be the case, otherwise the states are more 
or less bound. 
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and become independent of the specific profile. When line coin- 
cidences play a role, this might lead to certain errors though. 

4. The N in model ion 

We implemented a new Nitrogen model atom into the fastwind 
database, consisting of N n to N v. In the following, we provide 
some details of the Nm model ion, whilst the remaining ions 
will be described in Paper II. Our N m model consists of 41 lev- 
els, quite similar to the N in model as used within WM-basic. LS- 
coupled terms up to principal quantum number n — 6 and angu- 
lar momentum I = 4 have been considered. Table |B~T1 provides 
detailed information about the selected levels. All fine-structure 
sub-levels have been packed into one LS-coupled term. Q Two 
spin systems (doublet and quartet) are present and treated simul- 
taneously (see Fig. lB.ll for Grotrian diagrams). 

We account for all (193) allowed electric dipole radiative 
transitions between the 41 levels, as well as for (164) radia- 
tive intercombination transitions between the two spin system. 
Oscillator strengths have been taken from NIS10 when possible, 
and e lse from the WM -basic databaseQNIST N m data are mostly 
fr omlBell et al.l d!995 ). and from OPACITY project calculations 
bv lFernlev et al.l(ll999b . 

Roughly one thousand bound-bound collisional transitions 
between all levels are accounted for. (i) For all collisions be- 
tween the 1 1 lowest levels, 2s 2 2p, 2s 2p 2 , 2p 3 , and 2s 2 3/ (I - s, 
p, d), comprising doublet an d quartet terms, we u se the colli- 
sion strengths as calculat ed bvlStafford et al.l(l 19941) . from the ab 
initio R-matrix method (iBerrington et al.lll987l) . (ii) For most of 
the optically allowed transitions between higher levels (i .e., from 
level 1 1 as the lower one on), fhe lvan Regemorterld 19621) approx- 
imation is applied, (iii) For the optically forbidden transitions 
and the remainder of optically allowed ones (transitions in volv- 
ing th e highest level), the semi-empirical formula from lAllenl 
( 1 19731) (with Q = 1) is used. 

Photoionization cross-sections have been taken from 
the OPACITY Project on-line atomic database, TOPbaseQ 
(ICunto & Mendozal 19921). Th ese cross-sections have been com- 
puted bv lFernlev et al.l ([1999) by the R-matrix method using the 
close-coupling approximation and contain, as pointed out previ- 
ously, complex resonance structures. 

For excited N m levels with no OPACITY Project data avail- 
able (5g 2 G and 6g 2 G) and in those cases where we apply the 
explicit method for dielectronic recombination, resonance-free 
cross-sections are used, provided in terms of the lSeatonl (Tl958) 
approximation, 

a(v) = aoWivo/yY + (1 -y6)(v /v)' 5+I ], (4) 

with do the cross-section at threshold vo, and /3 and s fit parame- 
ters from the WM-basic atomic database. For most cross-sections, 
a reasonable consistency between these and the OPACITY 

4 When calculating the final synthetic profiles we use, when neces- 
sary, the un-packed levels by assuming that n,/g, - with occupation 
number iij and statistical weight g, - is similar for each sub-level within 
a packed level, due to strong collisional coupling. 

3 http://www.nist.gov/physlab/data/asd.cfm, firstly described in 
Kelleher et al. (1999) 

6 see [Pauldrach et al. 1 1994) . In brief, the atomic structur e code su- 
perstructure dEissner & N ussbaumerl H969 : lEi ssner 1991) has been 
used to calculate all bound state energies in LS and intermediate cou- 
pling as well as related atomic data, particularly oscillator strengths in- 
cluding those for stabilizing transitions. 

7 http://cdsweb.u-strasbg.fr/topbase/topbase.html 



Table 1. Model grid used by Mih alas & Hummed (1 19731) and in 
our test series. 



Model T 


3233 


3540 


3533 


3740 


3735 


4040 


4035 


4540 


5040 


r eff (kK) 
l°g£ 


32.5 
3.3 


35.0 
4.0 


35.0 
3.3 


37.5 
4.0 


37.5 
3.5 


40.0 
4.0 


40.0 
3.5 


45.0 
4.0 


50.0 
4.0 


M 


1.0 


0.35 


1.82 


0.58 


3.16 


0.93 


5.3 


8.05 


17.5 



Notes. The mass-loss rates provided in the last row (in units of 
10~ s M Q yr~') refer to our tests of wind effects (Sect. [53l l alone. All 
other tests have been performed with negligible M. 



project data is found, if one compares the resonance-free con- 
tribution only, see Fig. [2] 

The most important dielectronic recombination and reverse 
ionization processes are treated by the explicit method (in par- 
ticular, recombination to the strategic 3d level). Corresponding 
atomic data (wavelengths and oscillator strengths of the sta- 
bilizing transitions) are from the WM-basic atomic database as 
well. Finally, the cr oss-sec t ions f or collisional ionization are de- 
rived following the lSeatonl d 19621) formula, with threshold cross- 
sections from WM-basic. 



5. Nm (emission) line formation 

In the following, we describe the results of an extensive test se- 
ries regarding our newly developed N in model ion. In particu- 
lar, we check if the triplet appears in emission in the observed 
parameter range, and if the A4097 line remains always in ab- 
sorption. Note that these lines should be strongly correlated, i.e., 
the stronger the emission in the triplet, the weaker the absorp- 
tion at /14097, since both transitions share the level 3p (Fig. [TJ>. 
A change in the corresponding departure coefficient leads to a 
change in both lines. For example, if bj, p becomes diminished 
due to a more efficient drain by the 'two electron transitions' (see 
Sect. O, this leads to more emission at /U4634 - 4640 - 4642 
and to less absorption at /14097, due to less cascading. 

For all tests, we calculated model-grids that cover the same 
stellar parameters (O-type dwarfs and (super-) giants) as used by 
MH, listed in Table Q] All tests have been performed by means 
of fastwind, using our complete Nitrogen model atom. 

5. 1 . Comparison with the results from MH 

First, we test if we are able to reproduce the MH-results for 
0((f)) and O(f) stars. To this end, we need to invoke (almost) 
identical conditions, regarding both atmospheric and atomic 
models. Thus, we modified our N m ionic model, replacing part 
of our new data with those used by MH ('mixed' ionic model). 
In particular, we replaced data for dielectronic recombination (to 
the three draining levels, 2p 2 ( 2 D, 2 S, 2 P), to level 3d, and to few 
higher important levels - #14, 17, 20 and 21, see Table lB.fT i. the 
oscillator strengths for the 'two electron transitions' (Table 0, 
and the photoionization cross-sections for all levels below 3d 
(the cross-sections for the latter did already agree). 

Consistent with the MH-models, a Nitrogen abundance of 
[N]0 = 8.18 was adopted. This is a factor of 2.5 larger than the 



8 [A] = log A/H + 12, with A/H the number density of element A 
with respect to Hydrogen 
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Fig. 3. Comparison of equivalent widths (EW) for N m /14097 and /M4640/42. All atmospheric models calculated with 'pseudo line- 
blanketing' (see text). Black diamonds: results from MH; red triangles: new fastwind calculations with 'mixed' N m ionic model 
(new level structure, but important transition data from MH); purple squares: new fastwind calculations with new Nm model; blue 
asterisks: as squares, but DR-contribution to 3d level diminished by a factor of two. Here and in the following, positive and negative 
EWs refer to absorption and emission lines, respectively. 



Table 2. Oscil lator strengths for the 'two electron transition' as 
used by Mihalas & Hu mmer! ( 1 1973b and within our new atomic 
model. 



Transition 


fin (MH) 


fi u (used in this work) 


2p 2 2 D - 3p 2 P U 


2.50 ■ 10- 4 


4.38 ■ 10- J 


2p 2 2 S - 3p 2 P° 


2.00 ■ 10~ 4 


1.31 • 10~ 2 


2p 2 2 P - 3p 2 P° 


2.40 ■ 10- 2 


4.02 • 10~ 4 



solar one, [N Q ] = 7.78 (lAsplund et al .1120 05 )FI To account for the 
absence of a wind in their models, a negligible mass-loss rate of 
was used. 

Though line-blocking/blanketing could not be included into 
the atmospheric models in 1973, MH realized its significance 
and tried to incorporate some important aspect by means of a 
'pseudo-blanketing' treatment. The Nm photoionization edge 
(at 261 A) lies very close to the Hen ground state edge (at 
228 A) in a region of low continuum opacity and high emergent 
flux (if no line-blocking due to the numerous EUV metal-lines 
is present), and would cause severe underpopulation of the N m 
ground state if the blanketing effect is neglected, as shown by 
BM7 1 . MH argued that heavy line-blanketing and much lower 
fluxes are to be expected in this region. To simulate these effects, 
they extrapolated the He n ground-state photoionization cross- 
section beyond the Lyman-edge up to the N m edge, using a v~ 3 
extrapolation. To be consistent with their approach, we proceed 
in the same way by including this treatment into fastwind (see 
Fig.©. 



9 lAsplund et alj j2009t) provide a slightly larger value, [N Q ] = 
7.83 ± 0.05, where this difference is irrelevant in the following context. 
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Fig. 4. Comparison of fastwind Nm line profiles from model 
'T3740' with 'pseudo line-blanketing', using different atomic 
data sets. Solid (black): 'mixed' N m model (see text); dotted 
(red): new Nm model; dashed (blue): new Nm model, but DR- 
contribution to 3d level diminished by a factor of two. 

Figure [3] displays the comparison between the resulting 
equivalent widths from the MH models (black symbols) and our 
models using the conditions as outlined above (red triangles). 
Overall, the agreement for /U4634-4640-4642 is satisfactory^ 
and slight differences are present only for the hottest models. In 
agreement with the MH results, our profiles turn from absorp- 
tion into emission around T e ff ~ 37, 000 K for dwarfs and at 



10 Here and in the following, we only display the total equivalent 
widths of the /1/14640/4642 components; the behavior of /14634 is anal- 
ogous. 
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Fig. 6. EUV radiation temperatures for model 'T3735' with 
pseudo (dashed) and full (solid) line-blanketing. Important ion- 
ization edges are indicated as vertical markers. From left to right: 
N iv, He ii, N in (ground-state), and N in 2p 2 2 D. 



r eff ~ 33, 000 K for (super-)giants. As well, the ,14097 line is al- 
ways in absorption throughout the grid, though our calculations 
predict moderately more absorption in this line. All lines show 
the same trend in both sets, and the remaining differences might 
be attributed to still somewhat different atomic dataF*l 

We note already here that in all cases the emission is 
more pronounced in low-gravity objects. In high-gravity objects 
(dwarfs) part of the emission is suppressed because of higher 
collisional rates (oc « e ), driving the relative populations towards 
LTE. 

After demonstrating that we can (almost) reproduce the pro- 
files calculated by MH when similar conditions are applied, the 
next step is to investigate the effect of the new N m atomic data 
implemented during the present work. In fact, this leads to much 
more triplet emission, see Figs. [3] (purple squares) and [4] (dot- 
ted profiles). Even for the coolest models, where MH still ob- 
tain weak absorption, our calculations result in strong emission. 
This big difference is produced by larger DR rates into 3d and a 
larger drain of 3p by the 'two electron transitions', due to larger 
oscillator strengths (see Table ISliF 2 ! Note also that the /14097 line 
becomes weaker, in accordance with the correlation predicted 
above. 

The reaction of the emission strength on the dielectronic data 
also allows us to check the validity of the implementation of di- 
electronic recombination by the explicit method, described in 
Sect. IA.1I particularly the dependence on the oscillator strength 
of the stabilizing transition(s). Figures [3] (blue asterisks) and [4] 
(dashed) show the reaction of the triplet lines if we diminish the 
oscillator strength of the strongest stabilizing transition to 3d by 
a factor of two. This reduction leads to a significantly weaker 
emission throughout the whole grid, and a slight increase in the 
absorption strength of /14097. For a further test on the consis- 
tency between implicit and explicit method, see Appendix IA. 31 



11 number of levels, collisional data and LTE assumption concerning 
the quartet system levels by MH. 

12 A similar effect has already been noted by MH when they increased 
the corresponding oscillator strengths in their atomic model. 



5.2. Models with full line-blanketing 

In their study, MH could not consider the problem for realis- 
tic atmospheres accounting for a consistent description of line- 
blocking/blanketing, simply because such atmospheric models 
did not exist at that time. To investigate the differential effect on 
the N m emission lines, we calculated the same grid of models, 
now including full line-blanketin^Qas incorporated to fastwind, 
and compare with the MH calculations (Fig. [5] purple squares vs. 
black diamonds, respectively). 

Astonishingly, the triplet emission almost vanishes through- 
out the whole grid. This dramatic result points to the importance 
of including a realistic treatment of line-blanketing when inves- 
tigating the emission line problem in Of stars. It further implies, 
of course, that the new mechanism preventing emission in line- 
blanketed models needs to be understood, and that an alternative 
explanation/modeling for the observed emission must be found. 

Let us first investigate the 'inhibition effect', by consid- 
ering model 'T3735' that displays one of the largest reac- 
tions. The most important consequence of the inclusion of line- 
blocking/blanketing is the decrease of the ionizing fluxes in the 
EUV. Figure[6]dispiays corresponding radiation temperatures as 
a function of wavelength, for the 'pseudo line -blanketed', sim- 
ple model constructed in analogy to MH (dashed) and the fully 
line-blanketed model (solid). 

Indeed, the presence of substantially different ionizing fluxes 
around and longward of the N m edge is the origin of the dif- 
ferent triplet emission, via two alternative routes. For part of 
the cooler models (not the one displayed here), significantly 
higher radiation temperatures in the N m continuum (A < 26 1 A) 
of the 'pseudo line-blanketed' models lead to a strong ground- 
state depopulation. Moreover, the ground-state and the 2p 2 levels 
are strongly collisionally coupled, from the deepest atmosphere 
to the line-forming region of the photospheric lines. Thus, the 
ground-state and the 2p 2 levels react in a coupled way. If the 
ground-state becomes depopulated because of a strong radiation 
field, the 2p 2 levels become underpopulated as well, giving rise 
to a strong drain from level 3p and thus a large source-function 
for the triplet lines. Corresponding models including a realistic 
line-blocking (with lower EUV-fluxes) cause much less depop- 
ulation of the ground-state and the 2p 2 levels, thus suppressing 
any efficient drain and preventing strong emission in the triplet 
lines. 

For hotter models, e.g., model 'T3735' as considered here, 
the operating effect is vice versa, but leads to the same result. 
In the simple, MH-like model, the 2p 2 levels become strongly 
depopulated in a direct way, because (i) the ionizing fluxes at 
the corresponding edges (A > 261 A) are extremely high (the 
'pseudo line-blanketing extends 'only' until 261 A), and (ii) 
because also a strong ionization via doubly excited levels is 
present, again due to a strong radiation field at the corresponding 
transition wavelengths around 290-4 15 A. Due to the depopula- 
tion of the 2p 2 levels, the coupled ground-state becomes depop- 
ulated as well. 

On the other hand, the depopulation via direct and 'indirect' 
ionization from 2p 2 does not work for consistently blocked mod- 
els, because of the much lower ionizing fluxes. Insofar, the ex- 
tension of the He n ground-state ionization cross section by MH, 
to mimic the presence of line-blocking, was not sufficient. They 
should have extrapolated this cross-section at least to the edge 
of the lowest 2p 2 level. In this case, however, they would have 

13 Using backg round metallicities corresp onding to the 'older' solar 
abundances from Grevesse & Sauvall (l998), but keeping [N] = 8.18. 
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Fig. 5. Comparison of equivalent widths (EW) for Nm /14097 and /i/14640/42. Diamonds and triangles as Fig. [3] Purple squares: 
new results using new atomic data and realistic line-blanketing; blue asterisks: as squares, but DR-contribution to 3d level set to 
zero. 



found much lower N m emission, too low to be consistent with 
observations. 

In Fig. [7] we compare the corresponding departure coeffi- 
cients. Black curves refer to the simple, 'pseudo line-blanketed' 
model, and red curves to the corresponding model with full line- 
blocking/-blanketing. Obviously, both models display strongly 
coupled ground and 2p 2 states in the formation region of photo- 
spheric lines (here: tr oss > 0.02), where these states are much 
more depopulated in the simple model. 

For the conditions discussed so far, a significant depopula- 
tion of 2p 2 is a prerequisite for obtaining strong emission lines 
in the optical: only in this case, an efficient drain 3p — > 2p 2 due to 
cascading processes can be produced. The consequence of dif- 
ferent draining efficiencies becomes obvious if we investigate 
the run of bj p (dashed). In the pseudo-blanketed model, this level 
becomes much more depopulated relative to 3d (dashed-dotted), 
leading to much stronger triplet emission than in the full blan- 
keted model. 

The impact of the different processes can be examined in de- 
tail if we investigate the corresponding net rates responsible for 
the population and depopulation of level 3p (Fig. [8] upper panel). 
In this and the following similar plots, we display the dominating 
individual net rates (i.e., nJlji - n-Jiij > for population, with 
index i the considered level) as a fraction of the total population 
rateEl 

Indeed, there is a dramatic difference in the net rates that de- 
populate level 3p via the 'two electron' drain (dashed). Whereas 
for the 'pseudo line-blanketed' model the net rate into 2p 2 is the 
dominating one (resulting from the strong depopulation of 2p 2 ), 
this rate almost vanishes for the model with full line-blanketing. 
In contrast, the other two important processes, the population 
by 3d (solid) and the depopulation into 3s (dashed-dotted) are 
rather similar. Consequently, 3p becomes less depopulated com- 
pared to 3d in the fully blanketed models. Moreover, since the 



14 which in statistical equilibrium is identical to the total depopulation 
rate. 



1.2 r 




r. LlliJ _ 1 _ 1 _, '. 

100.000 10.000 1.000 0.100 0.010 0.001 

T Ross 

Fig. 7. Departure coefficients for important levels regarding the 
triplet emission, for model 'T3735' with pseudo (black) and full 
(red) line-blanketing (level 2p is the ground-state). The forma- 
tion region of the triplet lines is indicated by vertical dashed 
lines. See text. 



(relative) depopulation into 3s remains unaffected, and 3p has a 
larger population (see Fig. |7), level 3s becomes stronger popu- 
lated as well. Thus, the blanketed models produce more absorp- 
tion in /14097 (cf. Fig. [5] with Fig. [3] squares in upper panels). 

Since in the older MH-models the triplet emission is also due 
to the strong population of 3d via DR (Fig. [8] lower panel, red 
solid), we have to check how the presence of full line-blocking 
affects this process. In this case, the net-rate from DR (green 
solid) becomes even negative in part of the line forming region, 
i.e., the ionization via intermediate doubly excited states partly 
outweighs the dielectronic recombination. This difference origi- 
nates from two effects: (i) the photospheric electron densities in 
the older models are larger, because of missing radiation pres- 
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Fig. 8. Fractional net rates to and from level 3p (upper panel) and 
level 3d (lower panel), for model 'T3735' with pseudo (red) and 
full (green) line-blanketing. The formation region of the triplet 
lines is indicated by vertical dashed lines. 



sure from metallic lines. Higher electron densities imply higher 
recombination rates, (ii) In the fully blanketed models, the radi- 
ation field at the frequency of the main stabilizing transition to 
3d (777 A) is slightly higher (see Fig. |5), which leads to more 
ionization. Consequently, DR plays only a minor (or even op- 
posite) role in the fully blocked models, contrasted to the MH 
case. In the blocked models, the major population is via the 
4f level (dashed-dotted). Note also that the Swings mechanism, 
i.e., a population via the resonance line at 374 A (dashed), does 
not play any role in photospheric regions, as already argued by 
BM71 and shown by MH. 

For a final check about the importance of DR to level 3d in 
models with full line-blanketing, we tested its impact by switch- 
ing off the main stabilizing transition. From Fig. [5] (blue aster- 
isks), it is obvious that the effect indeed is marginal throughout 
the whole grid. Thus, a correct treatment of line-blocking seems 
to suppress both the efficiency of the draining transitions and the 
dielectronic recombination, and we have to ask ourselves how 
the observed triplet emission is produced, since the presence of 
line-blocking cannot be argued away. 

Before tackling this problem, in Fig. [9] w e compare our 
results (red triangles ) with those from tlusty dHubenvl tl988: 
iHubenv & Land Il995lf^l (blue asterisks), as provi ded by the 
OSTA R2002 gridB (La nz & Hubenvll2003l: see also iHeap et all 
2006). Contrasted to all other similar plots, we used an abun- 
dance of [N] = 7.92, to be consistent with the tlusty grid. 
Obviously, our fastwind results for the Nm emission lines in 
dwarfs match exactly those from tlusty, whereas our results for 
supergiants display somewhat less emission (or more absorp- 
tion). On the other hand, the absorption line at A4Q91 is sys- 



15 a code that assumes plane-parallel geometry, hydrostatic and ra- 
diative equilibrium, and calculates line-blanketed NLTE model atmo- 
spheres and corresponding synthetic profiles. Due to its restrictions, 
only objects with negligible winds can be analyzed. 

16 Equivalent widths from own integration. 
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Fig. 10. Comparison of equivalent widths (EW) for N m /14097 and /L14640/42. Black diamonds as in Fig.[3](results from MH). Red 
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tematically stronger in all our models, which points to different 
oscillator strengths. To rule out potential differences in the at- 
mospheric stratification, we calculated an additional grid with 
temperature and electron structure taken from tlusty, smoothly 
connected to the wind structure as calculated by fastwind. 
Corresponding results are displayed by purple squares, and differ 
hardly from those based on the original fastwind structure, ex- 
cept for model T3233 where the differences in T and n e are larger 
than elsewhere. The remaining discrepancies for supergiants can 
be explained by somewhat higher EUV fluxes in the tlusty mod- 
els, favouring more N m triplet emission. Nevertheless, the dis- 
agreement between fastwind and tlusty is usually much weaker 
than between these two codes and the results by MH, which un- 
derpins our finding that the triplet emission becomes strongly 
suppressed in fully blanketed models. 

5.3. The impact of wind effects 

So far, one process has been neglected, namely the (general) 
presence of winds in OB-stars. Note that MH had not the re- 
sources to reproduce O-stars with truly "extended" atmospheres. 
Regarding their O(f) and 0((f)) objects, on the other hand, there 
was no need to consider wind effects, because the observed N in 
triplet emission could be simulated by accounting for DR and 
'two-electron' drain alone. Actually, MH pointed out that in Of- 
supergiants (with denser winds) the Swings mechanism could 
play a crucial role in the overpopulation of 3d: Velocity fields 
are able to shift the resonance lines into the continuum, allowing 
them to become locally more transparent, and deviations from 
detailed balance and strong pumping might occur. 

With the advent of new atmospheric codes, we are now 
able to investigate the general role of winds, and to explain 
how the emission lines are formed within such a scenario. To 
this end, we have re-calculated the model grid from Table [1] 
now including the presence of a wind. Prototypical values for 



terminal velocity and velocity-field exponent have been used, 
Vco = 2,000 kms~', (3 = 0.9, and mass-loss rates were in- 
ferred from the wind-momentum l uminosity relationship (WLR) 
provided by Repolust et al. (2004), which differentiates between 
supergiants and other luminosity classes. Clumping effects have 
been neglected, but will be discussed in Sect. [8] For a summary 
of the adopted mass-loss rates, see Table [1] 

With the inclusion of winds now, we almost recover the 
emission predicted in the original MH calculations (Fig. ITOb . 
Less emission is produced only for the two hottest models, 
where the inclusion of a wind has no effect. For all other models, 
however, the wind effect is large, both for the supergiants and for 
the dwarfs with a rather low wind-density. 

To understand the underlying mechanism, we inspect again 
the fractional net rates, now for the wind-model 'T3735' 
(Fig. [m >. Obviously, the wind induces a significant overpopu- 
lation of the 3d level via the ground state rather than the dielec- 
tronic recombination, just in the way as indicated above. Due to 
the velocity field induced Doppler-shifts, the resonance line(s) 
become desaturated, the rates are no longer in detailed balance, 
and considerable pumping occurs because of the still quite large 
radiation field at 374 A and the significant ground-state popula- 
tion (larger than in the MH-like models). Interestingly, this does 
not only happen in supergiants, as speculated by MH, but also 
in dwarfs (at least those with non-negligible mass-loss rate), be- 
cause the velocity field sets in just in those regions where the 
emission lines are formed^ 

The impact of the wind on the population of the 3p level is 
not as extreme, and the 'two electron transition' drain remains as 
weak as for the wind-free model. Thus, the presence of emission 
relies mostly on the overpopulation of the 3d level. 



17 An analogous desaturation through velocity field ind uced Doppler 
shifts in stars with low mass loss rates was found by iNaiarro et al.l 
i ll 9961) in the He i resonance lines of B giants. 
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Fig. 11. Fractional net rates to and from level 3p (top) and level 
3d (bottom), for model 'T3735' (full line-blanketing) with no 
wind (red) and with wind (green). 



Again, we test the (remaining) influence of dielectronic re- 
combination by switching off the stabilizing transition to 3d. 
Though there is a certain effect, the change is not extreme. Note 
also that the net rate from dielectronic recombination is larger in 
the wind model than in the wind-free model, see Fig. QT| (green 
solid line). Nevertheless, we conclude that dielectronic recombi- 
nation plays, if at all, only a secondary role in the overpopulation 
of the 3d level when consistent atmospheric models are consid- 
ered. The crucial process is pumping by the resonance lines. 

The reaction of /14097 on velocity field effects is more com- 
plex. On the one side, there is still the cascade from 3p to 3s, 
giving rise to a certain correlation. On the other, the resonance 
line to 3s (at 452 A) becomes efficient now, and can either feed 
(as argued in Sect. [2]) or drain level 3s, in dependence of its op- 
tical depth. Under the conditions discussed here, the resonance 
line is pumping at cooler temperatures, with a zero net effect on 
the strength of ,14097 (since the cascade from 3p becomes some- 
what decreased, compared to wind-free models). At higher tem- 
peratures, the resonance line becomes optically thin{3 because 
of a lower ground-state population, and level 3s can cascade to 
the ground-state. Thus, the absorption strength of ,14097 might 
become significantly reduced, which explains, e.g., the strong 
deviation of model 'T4035' from the MH predictions. For this 



Table 3. Stellar and wind parameters of our model grid used 
to compare synthetic N n/N m profiles from fas twind and cmf- 
gen. T he grid is a subset of the grid presented bv lLenorzer et alj 
(2003). 



Luminosity class V 



Model 


Teff 


R, 


l°gg 


M 




P 




(K) 


(Re) 


(cgs) 


(lO^Moyr 1 ) 


(kms" 1 ) 




d2v 


46100 


11.4 


4.01 


2.52 


3140 


0.8 


d4v 


41010 


10.0 


4.01 


0.847 


2850 


0.8 


d6v 


35900 


8.8 


3.95 


0.210 


2570 


0.8 


d8v 


32000 


8.0 


3.90 


0.056 


2400 


0.8 


dlOv 


28000 


7.4 


3.87 


0.0122 


2210 


0.8 


Luminosity class I 


Model 


Teff 


R, 


l°gg 


M 




P 




(K) 


(Re) 


(cgs) 


(lO^Meyr" 1 ) 


(kms" 1 ) 




s2a 


44700 


19.6 


3.79 


12.0 


2620 


1.0 


s4a 


38700 


21.8 


3.57 


7.35 


2190 


1.0 


s6a 


32740 


24.6 


3.33 


3.10 


1810 


1.0 


s8a 


29760 


26.2 


3.21 


1.53 


1690 


1.0 


slOa 


23780 


30.5 


2.98 


3.90 


740 


1.0 



and the hotter models, ,14097 is very weak, and even appears in 
(very weak) emission for model 'T5040'. 



6. Comparison with results from cmfgen 

In this section, we compare the results from our fastwind models 
with corresponding ones from cmfgen, a code that is considered 
to produce highly reliable synthetic spectra, due to its approach 
of calculating all lines (within its atomic database) in the co- 
moving frame. For this purpose, we used a grid of models for 
dwarfs and supergiants in the O and early B-sta r range, which 
has a l ready been used in previous compariso ns dLeno rzer et al. 
120041: iPuls eTal] 120051: iRepolust et alj|2005l) . fastwind models 
have been calculated with three 'explicit' atoms, H, He and our 
new N atom. Stellar and wind parameters of the grid models 
are listed in Table [3] with wind parameters following roughly 
the WLR for Galactic stars. The model designations correspond 
only coarsel y to spectral types and are, with respect t o recent 
calibrations dRepolust et al.i r2004: Ma rtins et alJl2005a) . some- 
what too early. All calculations have been performed with the 
'old' solar Nitrogen abundance, [N] = 7.92 (Greves se & Sauval 
Il998h . and a micro-turbulence v tul -b = 15 kms . 

For the comparison between fastwind and cmfgen results, we 
consider useful diagnostic N n (to check the cooler models) and 
N m lines in the blue part of the visual spectrum. Table|4]lists im- 
portant N in lines in the range betwee n 4,000 to 6,500 A, together 
with multiplet numbers according to iMoord d 19751) . to provide 
an impression of how many independent lines are present and 
which lines belong to the same multiplet[3 Additionally, we pro- 
vide information about adjacent lines present in the spectra of 
(early) B and O stars, where the major source of contamination 
arises from O n and C m. 

This set of lines can be split into five different groups. Lines 
that belong to the first group (#2,5,9-1 1 in Tabled are produced 
by cascade processes of the (doublet) series Is 2 2s 2 nl with n = 
3,4,5 and / = s, p, d, f, g, and additional (over-) population of 



18 Note that the oscillator strength is more than a factor of 10 lower 
than that of the resonance line feeding 3d. 



19 A similar table regarding N n, N iv and N v lines will be presented 
in Paper II. 
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the 3d level. As outlined in the previous sections, the behavior 
of these lines is strongly coupled. 

The second group (lines #6-8), AA45 10 -45 14 -45 18, results 
from transitions within the quartet system. We consider only the 
three strongest components of the corresponding multiplet, that 
are also the least blended ones. 

The third group (lines #1,3,4) represents three lines, at 
,1,14003, 4195, and 4200 A, that are formed between higher lying 
levels within the doublet system. These lines are weaker than the 
ones from the previous two groups, but still worth to use them 
within a comparison of codes and also within a final abundance 
analysis. Note that N mvl4195 and N ni/14200 are located within 
the Stark- wing and the core of He n/14200, respectively, which 
requires a consistent analysis of the total line complex. 

Lines at ,1,15320, 5327 A (#12-13) and /L16445, 6450, 6454, 
6467 A (#14-17) comprise the fourth and fifth group, respec- 
tively. The former set of lines is located in a spectral region that 
is rarely observed, and the latter comprises a multiplet from the 
quartet system, in the red part of the visual spectrum. 

Although most of the lines listed in Table |4] are (usually) 
visible in not too early O-type spectra, their diagnostic po- 
tential for abundance determinations is different. The lines at 
/L14510-4514-4518 (from the quartet system) are certainly the 
best candidates to infer abundances, since they are quite strong 
and their formation is rather simple. Also the N m triplet itself 
provides valuable information. Due to its complex formation 
- when in emission - and additional problems (see below and 
Sect [7]l, these lines should be used only as secondary diagnos- 
tics, whenever possible. The remaining lines in the blue part of 
the visual spectrum are rather weak and/or strongly blended with 
adjacent lines (see Table©. In particular, the (theoretically) very 
interesting transition 3p — > 3s, Nni/14097, is located within the 
Stark- wing of Hg. Thus, these lines should be used preferentially 
as a consistency check, and employed as a direct abundance indi- 
cator only at low rotation. Lines in the yellow part of the optical 
(group four) have not been considered by us so far, since our 
observational material does not cover this spectral range, and 
we are not able to judge their diagnostic value for O-type stars. 
Finally, lines from group five in the red are usually rather weak, 
and might be used only in high S/N spectra of slowly rotating 
stars. 

A detailed comparison between the various N n (for the later 
subtypes) and N m lines from fastwind and cmfgen is provided 
in Appendix[C] Overall, the following trends and problems have 
been identified. 

For models dl(8v, d8v, and slOa (Figs.|Cj]E2] and|C3] re- 
spectively), the agreement of the N n lines is almost perfect. For 
model s8a (Fig. IC.4D . on the other hand, big discrepancies are 
found. Most of our lines are much stronger than those from cm- 
fgen, because of the following reason. Within the line-forming 
region, the electron density as calculated by fastwfnd is a factor 
of = 8 higher than the one as calculated by cmfgen, thus enforc- 
ing higher recombination rates from N in to N n, more N n and 
thus stronger lines. This is the only model with such a large dis- 
crepancy in the electron density (e.g., the electron densities from 
models slQa and s6a agree very well), and the reason for this 
discrepancy needs to be identified in future work. Nevertheless, 
this discrepancy would not lead to erroneous Nitrogen abun- 
dances: When analyzing the observations, a prime diagnostic 
tool are the wings of the Hydrogen Balmer lines that react, via 
Stark-broadening, sensitively on the electron density. After a fit 
of these wings has been obtained, one can be sure that the elec- 
tron stratification of the model is reasonable (though the derived 



Table 4. Diagnostic N m lines in the optical, together with adja- 
cent blends. 



# Wavelength M# low-up Blends 

(A) 

1 4003.58 17 20-33 On ,14007.46 

2 4097.33 1 8-10 On ,14097.26, 4098.24, H, ,14101.74 

3 4195.76 6 13-22 On ,14192.52, 4196.26, Si m,14195.59, 
Hen ,14200.00 

4 4200.07 6 13-22 Hen ,14200.00 

5 4379.11 18 21-34 On ,14378.03, 4378.43, Cna4379.47, 
Nn ,14379.59 

6 4510.88 3 12-16 N in ,14510.92, Nen,t4511.42 

7 4514.86 3 12-16 Om ,14513.83, Ne il14514.88, 
Cm,14515.81, 4516.77 

8 4518.14 3 12-16 Nen ,14518.14, Om,14519.62 

9 4634.14 2 10-11 Siiv ,14631.24, Orv,14632 

10 4640.64 2 10-11 On ,14638.86, Si nU4638.28 

11 4641.85 2 10-11 On ,14641.81, 4643.39, Nn,14643.08 

12 5320.82 21 22-30 On ,15322.53 

13 5327.18 21 22-30 - 

14 6445.34 14 19-25 - 

15 6450.79 14 19-25 Civ ,16449.90 

16 6454.08 14 19-25 On ,16457.05, Nn,16457.68 

17 6467.02 14 19-25 - 

Notes. Line numbers refer to the transi tions indicated in Fig. IB. II 'M#' 
are the multiplet numbers according to M oord d 19751) . and 'low-up' re- 
fer to the corresponding (packed) lower and upper levels as provided in 
Table lO 



gravity might be erroneous then). With a 'correct' electron den- 
sity, however, the abundance determination via N n should be no 
further problem. 

Let us concentrate now on the N m lines, beginning with the 
triplet ,1,14634 - 4640 - 4642 and accounting for the fact that 
at cooler temperatures the ,14642 component is dominated by an 
overlapping On line (see Table [4]). Except for the dwarf model 
d6v and the supergiant models s8a and s6a, both codes pre- 
dict very similar absorption (for the cooler models) and emission 
lines, thus underpinning the results from our various tests per- 
formed in Sect. [5] Interestingly, our emission lines for the hottest 
dwarfs (d4v and d2v, Figs. IC.8I and IC.9b are slightly stronger 
than those predicted by cmfgen, a fact that is not too worrisome 
accounting for the subtleties involved in the formation process. 

What is worrisome, however, is the deviation for models 
d6v, s8a and s6a (Figs. |C7] ICTTl and ICL21 . Whereas fast- 
wind predicts only slightly refilled absorption profiles for all 
three models, cmfgen predicts almost completely refilled (i.e., 
EW « 0) profiles for d6v and s8a, and well developed emission 
for s6a. 

We investigated a number of possible reasons for this dis- 
crepancy. At first, the difference in electron density for model 
s8a (see above) could not explain the strong deviations. At least 
for our fastwind models, dielectronic recombination does not 
play any role, and we checked that this is also the case for the 
cmfgen models, by switching off DR. As it finally turned out, the 
physical driver that overpopulates the 3d level in cmfgen in the 
considered parameter range is a coupling with the O in resonance 
line at 374 A, which is discussed in the next section. 

With respect to N m,i4097, we find a certain trend in the de- 
viation between fastwind and cmfgen absorption profiles. Both 
for dwarfs and supergiants, our line is weaker at cooler tempera- 
tures, similar at intermediate ones and stronger at hotter temper- 
atures. 



4003.58 


17 


20-33 


4097.33 


1 


8-10 


4195.76 


6 


13-22 


4200.07 


6 


13-22 


4379.11 


18 


21-34 


4510.88 


3 


12-16 


4514.86 


3 


12-16 


4518.14 


3 


12-16 


4634.14 


2 


10-11 


4640.64 


2 


10-11 


4641.85 


2 


10-11 


5320.82 


21 


22-30 


5327.18 


21 


22-30 


6445.34 


14 


19-25 


6450.79 


14 


19-25 


6454.08 


14 


19-25 


6467.02 


14 


19-25 



O ii ,14007.46 

On ,14097.26, 4098.24, H a ,14101.74 
On ,14192.52, 4196.26, Sim,14195.59, 
He ii ,14200.00 
He ii ,14200.00 

On ,14378.03, 4378.43, C m,14379.47, 
Nn ,14379.59 
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For N ni/14379, the cmfgen models predict more absorption 
for models with r e ff < 35, 000 K, because of an overlapping O n 
and C m line. For hotter models, the presence of O n is no longer 
important, but there is a still a difference for /14379, going into 
emission at the hottest temperatures. 

For the quartet multiplet at /U4510 - 4514 - 4518, we find 
a similar trend as for /14097. For dwarfs, the lines are weaker at 
cooler temperatures, similar at d6v and stronger in absorption 
at d4v. Likewise then, the emission at d2v is weaker than in 
cmfgen. For supergiants, the situation is analogous, at least for 
all models but the hottest one where both codes predict identical 
emission. 

Finally, lines /14003 and /14195 (note the blueward Si m blend 
at hotter temperatures) show quite a good agreement for both 
dwarfs and supergiants, with somewhat larger discrepancies only 
for models d4v and s4a. 

In summary, the overall agreement between N n and N m pro- 
files as predicted by fastwind and cmfgen is satisfactory, and 
for most lines and models the differences are not cumbersome. 
Because of the involved systematics, however, abundance anal- 
yses might become slightly biased as a function of temperature, 
if performed either via fastwind or via cmfgen. Additionally, we 
have identified also some strong deviations, namely with respect 
to N ii line-strengths for supergiants around T e gm 30 kKriand 
regarding the emission strengths of the N m triplet, for models 
d6v, s8a, and s6a. 

7. Coupling with Om 

After numerous tests we were finally able to identify the origin 
of the latter discrepancy. It is the overlap of two resonance lines 
from O m and N m around 374 A (for details, see Table [5]l that 
is responsible for the stronger emission in cmfgen models com- 
pared to the fastwind results. Whereas this process is accounted 
for in cmfgen, our present fastwind models cannot do so, be- 
cause Oxygen is treated only as a background element, and no 
exact line transfer (including the ov e rlap) i s performed. We note 
that this is not exactly the iBowenl (Il935l) mechanism as men- 
tioned in Sect. [2] since this mechanism involves also the overlap 
between another O m resonance line and the He n Ly-a line at 
303 A, which does not play a strong role in our cmfgen mod- 
els as we have convinced ourselves^j Here, only the coupling 
between N m and O ra at 374 A is decisive. 

In the following we discuss some details of the process, by 
means of the cmfgen model s6a that displays the largest differ- 
ence to the fastwind predictions. For a further investigation, we 
'decoupled' Nm from Om by setting the corresponding oscil- 
lator strength of the transition Om 2p 2 3 P2 — > 3s 3 P ( ] ) to a very 
low value, and compared the results with those from the 'cou- 
pled' standard model. Figure [12] shows that the well developed 
emission lines from the standard models switch into absorption, 
and become very similar to the corresponding fastwind profiles 
(Fig. IC.12l i. On the other hand, discarding DR by neglecting the 
resonances in the photo cross-sections from level 3d had almost 
no effect on the profiles, in accordance with our previous argu- 
ments. (Actually, without resonances there is even more emis- 
sion than before, which shows that for this model the ionization 
via doubly excited levels outweighs DR.) 

Figure [13] displays the corresponding NLTE departure coef- 
ficients for the N in ground state, level 3p and level 3d. Because 

20 rooted within a different stratification of electron density. 

21 Insofar, the Om lines at X13340, 3444, 3759 mentioned in Sect. [2] 
remain in absorption. 



Table 5. Overlapping N m and O in resonance lines around 374 
A. 



transition(s) 


A 

(A) 


low 


gi 


up 




*/ 


ratio 


Nm (packed) 




2p *P<> 


6 


3d 2 D 


10 


2.5 




N m coupl. comp. 


374.434 


V^ 2 


4 


3d 2 D 5/2 


6 


1.5 


0.60 


O m (packed) 




2p z 3 P 


9 


3s J P" 


9 


0.72 




O m coupl. comp. 


374.432 


2p 2 3 P 2 


5 


3s 3 P° 


3 


0.06 


0.14 



Notes. Packed levels/transitions and overlapping ('coupled') ones, 
'low' = lower term/level, 'up' = upper term/level, g\ and g u correspond- 
ing statistical weights. 'Ratio' indicates the opacity ratio between over- 
lapping and packed lines. Individual components at 374.20. . . 374.44 A 
for Nm (3 components) and at 373.80. .. 374.43 A for Om (6 com- 
ponents). Note that 0.01 A correspond to 8kms~'. Data from NIST 
database. 
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Fig. 12. N m AA4634 - 4640 - 4642 for three versions of model 
s6a, as calculated by cmfgen. Solid black: standard version; 
Dotted/red: photo cross-sections for level 3d without resonances 
(i.e., no DR); dashed/blue: N in resonance line at 374 A forbid- 
den to be affected by O m. For the latter model, the optical triplet 
lines remain in absorption! 

of the superlevel approach, there is no differencial between level 
3d 2 Ds/2 (coupled to O in) and 3d 2 Di/2 (not coupled). Whilst the 
ground state and the 3p level remain similar for all three models 
(standard, 'decoupled' and 'coupled' with resonance-free photo 
cross-section for 3d), level 3d becomes much stronger overpop- 
ulated in the two 'coupled' models, compared to the 'decoupled' 
one. 

The origin of this stronger overpopulation becomes obvious 
from Fig. [14] The upper panel displays the source functions of 
the N m/O in resonance lines for the 'decoupled' and 'coupled' 
case (again, no difference for individual components, due to the 
superlevel approach). For the 'decoupled' case, the N m source 
function is significantly lower than the O m one, whereas for the 
'coupled' standard model both source functions are identical, at 
a level very close to the 'decoupled' O m case. These changes 
are visualized in the lower panel: because of the line overlap, 
the Nm source function increases by a factor up to 1.35 (in the 



assumed to be balanced by collisions. 
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Fig. 13. NLTE departure coefficients for three versions of model 
s6a, as calculated by cmfgen. Lower group of curves: Nm 
ground state; intermediate group: level 3p, upper group: level 
3d. See text. 

line forming region of the optical triplet lines), whereas O in de- 
creases only marginally. 

The source function equality itself results from the strong 
coupling of the two resonance lines and the fact that both of 
them are optically thick from the lower wind on, with similar 
opacities (proportional to the product of cross-section, ionization 
fraction and abundance) It can be shown that the magnitude of 
the coupled source function mostly depends on those processes 
that determine the individual source functions in addition to the 
mean line intensity. This can be easily seen by using the Sobolev 
approximation, 



J = (l-p)S +/3J C , 



(5) 



with scattering integral J, loc al and core es cape probabilities p 
and yS c , and core intensity I c (ISobolevlll960l) . In case of over- 
lapping lines, the source function S and the opacities need to 
account for all components with appropriate weights. 

We now assume that the individual source functions (for the 
packed levels) can be approximated by 



J + Sq, 



(6) 



where 6^ and So correspond to the additional source terms 
(mostly from cascades to the upper levels), and the well-known 
factofl in front of J has been approximated by unity. 

J includes the contributions from all transitions between the 
fine structure components, 7=2 gfiJil X! gfi, exploiting the fact 
that the reduced occupation numbers for the fine structure levels 
are equal, m/gi = nj/gj. 

Let us first consider the case of no fine structure splitting, i.e., 
the overlapping resonance lines should be the only ones connect- 
ing the upper and lower levels. Without line overlap, we obtain 
the well-known result (solving for Eqs.|5]and|6]in parallel) 



, decoup _ A:,N^c + <?N 



Recoup _ A:,Q^c + #0 



N 



Pn 



'O 



Po 



(7) 



23 In the formation region of the optical lines, the opacity ratio be- 
tween the N in and O m resonance lines ranges from roughly 5 for the 
coolest models to 0.5 for the hottest ones. 

24 corresponding to (1-e) in a two-level atom, with s roughly the ratio 
of collisional to spontaneous radiative rate coefficient, for the downward 
transition. 




Fig. 14. Line source functions for the N m/O m resonance lines at 
374 A, for two versions of model s6a, as calculated by cmfgen. 
Upper panel: Nm 'decoupled' (dashed/blue); Om 'decoupled' 
(dashed-dotted/red); Nm and Om 'coupled' (solid, identical). 
Lower panel: Source function ratios. Solid: N m ('coupled')/Nm 
('decoupled'); dashed: Om ('coupled')/Om ('decoupled'). The 
optical triplet lines are formed in the range indicated by the dot- 
ted lines. 



Likewise, but accounting for the line overlap and assuming 
P <K 1 for both components, we find 



coup 
N 



c, coup 



PJc 
P 



(J N <5o 
+ h 



r, decoup 



Pn Po 

decoup 



(8) 



Thus, the coupled source functions are larger than in the decou- 
pled case, but (almost) independent on the ratio of their contri- 
bution. (In the above equation, the first term of the sum does not 
depend on any specific opacity as long as the lines - coupled or 
decoupled - remain optically thick). 

For the case of two multiplet lines from both N m and O m, 
with one of each overlapping, the corresponding result for the 
packed levels reads 



t coup 



^.coup 



p + 3\8 N + Po'' 



(9) 
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where the escape probabilities /?n and f3o refer to the total opac- 
ity arising from both multiplet lines, and the factor '1/3' can be 
traced down to the fact that three (optically thick) lines partici- 
pate in total, two single lines from N m and O m, and one coupled 
line complex. 

The corresponding source functions for the decoupled case 
with packed levels would read 

r,decoup _ ficlc 1 C>N „decoup _ fic^c 1 <^0 fir\\ 

5 n " + 2A? S ° = T Wo' (10) 

(the factor of two arising from two participating multiplet lines), 
and a comparison with Eq. [9] shows that in most cases the cou- 
pled source function would lie in between the corresponding 
ones that neglect the line overlap. Again, our result is (almost) 
independent of the opacity ratio between the overlapping N m 
and O in lines but also independent of the weights of the individ- 
ual lines within each multipletPI 

Generalization to more multiplet lines is straightforward, and 
our analytic result compares well to the actual case where the 
N in and O in source functions for the coupled case are identical 
and lie in between the values for the decoupled situation, see 
Fig- [HI upper panel. 

Thus, whenever the source function of the O iii/1374 line is 
significantly larger than the one from N in, a decent effect on the 
emission strength of the optical triplet lines is to be expected. 
This situation is particularly met around T e ff = 30 to 33 kK, since 
in this region the upper level of the O in line is significantly pop- 
ulated by cascades from higher levels. 

First test calculations with a simplified Oxygen model atom 
performed by fastwind confirm the general effect, but realistic 
results cannot be provided before a detailed model atom has been 
constructed. Let us note, however, that most of the other optical 
N in lines are barely affected by the resonance line coupling, and 
that these lines can be used for diagnostic purposes already now. 

In Paper II we derive Nitrogen abundances for LMC O-stars 
from the VLT-FLAMES survey of massive stars. Though there 
are only few objects in the critical temperature range, at least 
for one object, Nl 1-029 (09.71b), we encountered the problem 
that the observed, refilled Nm triplet (EW » 0) could not be 
reproduced by fastwind, though with cmfgen. We interpret this 
problem as due to the resonance line overlap, but stress also the 
fact that other diagnostic lines enable a satisfactory abundance 
analysis. 

8. Influence of various parameters 

Let us finally investigate the influence of important parameters 
on the strength of the optical emission lines. We stress that the 
following results have an only qualitative, differential character, 
as long as the coupling with O in has not been accounted for, at 
least in the range T e fi S 35kK. 

Nitrogen abundance. Figure[15]displays the reaction on a vari- 
ation of Nitrogen abundance and mass-loss rate. All models 
have been calculate d with background abundances following 
Asplund et al. (2005). The (purple) squares correspond (almost, 
except for the background) to our previous results for [N] = 
8.18 (0.4 dex larger than the solar value), and mass-loss rates 
according to TableQ] Reducing the Nitrogen abundance to solar 
values, [N] = 7.78, results in considerably less emission (black 
diamonds), which is a consequence of the fact that the relative 

25 as long as /? <sc 1. 



overpopulation of the 3d level decreases significantly when the 
formation depth proceeds inwards. The influence on N ni/14097, 
on the other hand, is less extreme, and follows the standard trend 
that a higher abundance results in stronger absorption lines. 

Mass-loss rates. Also in Fig. [15] asterisks (for [N] = 8.18) and 
triangles (for [N] = 7.78) correspond to a situation where the 
mass-loss rates have been decreased, by a factor of two com- 
pared to Table [1] Whereas the effect for the enriched models is 
significant, the models with solar nitrogen abundance are much 
less affected, though the predicted emission is still much stronger 
than for models without a wind at all, cf. Fig.fTOl(triangles). 

The origin of this different behaviour is, again, rooted in 
the specific run of the relative overpopulation of level 3d. For 
a large Nitrogen abundance, the formation takes place in a re- 
gion where the overpopulation increases quickly with increasing 
mass-loss rate, whereas for a lower abundance (deeper forma- 
tion) it remains rather constant over a certain range of formation 
depths, such that a moderate change in M has much less effect. 
Figures Q~3] and [14] (upper panel), with a formation region cor- 
responding to [N] = 7.92, visualize the general situation (irre- 
spective of whether coupling with O m is important or not): an 
increase of the Nitrogen abundance shifts the formation depth 
to the right (towards lower tr oss ), where the departure coeffi- 
cient of level 3d and the source function display a 'bump', due 
to the large velocity gradient, inducing very strong pumping by 
the resonance line(s). A decrease of the wind-strength, on the 
other hand, shifts the position of the bump to the right, while 
(almost) preserving the formation depth with respect to tr oss . 
Thus, in case of a large Nitrogen abundance the strong overpop- 
ulation is lost when M decreases, leading to a corresponding loss 
of emission strength. For a lower abundance, however, the for- 
mation takes place in a region where the overpopulation displays 
a 'plateau', and a reduction of M, i.e., a shift of the bump to the 
right, has a much weaker effect. 

Wind clumping. Since the Nm emission lines of O-stars are 
formed in the middle or outer photosphere|3 wind clump- 
ing has no direct effect on their strength, as long as the pho- 
tosph ere remains unclumped (which is most likely the case , 
e.g. JPuls et alJl2006HSundqvist et aD201 lLlNaiarro et al.ll201 ll) . 
Only if the transition region wind/photosphere were clumped, 
such direct effects are to be expected, due to increased recom- 
bination and optically thick clumps affecting the resonance lines 
dSundqvist et al.ll201 lh . Indirect effects, however, can be large, 
since clumpy winds have a lower mass-loss rate than corre- 
sponding smooth ones. Since typical M reductions are expected 
to be on the ord er of 2. . . 3 when acc ounting for micro- and 
macro-clumping ( Sundq vist et al.ll201ll and references therein), 
Figure[l5](asterisks vs. squares and triangles vs. diamonds) gives 
also an impression on the expected effects when clumping were 
in- or excluded in the spectrum synthesis. 

Background abundances. Figure [16] displays the effects when 
only the background abundances are varied, while, inconsis- 
tently, the Nitrogen abundance is kept at its solar value and the 
wind-strengths correspond to the Galactic WLR. Diamonds refer 
to solar composition, triangles to Z/Z Q = 0.5 (roughly LMC) and 
squares to Z/Z = 0.2 (roughly SMC). In line with our reason- 



26 only for a very large Nitrogen abundance or mass-loss rate, ('slash 
stars' or Wolf-Rayets) the formation takes place in the wind. 
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Fig. 15. Comparison of equivalent widths (EW) for Nm /14097 and /L14640/42 for models with different Nitrogen abundances 
and wind-strengths. Black diamonds: [N] = 7.78, mass-loss rates as in Table Q] Red triangles: as black diamonds, but with half 
the mass-loss rate. Purple squares: [N] = 8.18, mass -loss rates as in Ta ble Q] Blue asterisks: as purple squares, but with half the 
mass-loss rate. All background abundances are solar (Asplundetal. 200l). 
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Fig. 16. Comparison of equivalent widths (EW) for N m ,14097 and /1/14640/42 for models with different background abundances, 
[N] = 7.78 and mass-loss rates as in TableQ] Black diamonds: Z/Z - 1 .0; red triang les: Z/Z Q = 0.5; purple squares: Z/Z Q = 0.2. 
Background abundances scaled to solar abundance pattern from As plund et alJ d2005l) 



ing from Sect. 15.21 the emission increases significantly when the 
background abundances decrease, due to reduced line-blocking. 
Note that for objects with unprocessed Nitrogen this effect might 
be compensated by a lower Nitrogen abundance, a nd, even more, 
becau se of lower wind-strengths (M oc (Z/Z ) a7 , iMokiem et al.l 
120071) . Note also that the correlation between triplet emission 
and /14097 absorption strength has become rather weak, due to 
the influence of the resonance line connected to level 3 s (see 
Sect. I5.31 l. Let us finally stress that because of the strong de- 



pendence of emission line strength on line-blocking, a sensible 
assumption on the atmospheric iron abundance is required, due 
to its dominating effect on the EUV fluxes. 

9. Summary and conclusions 

In this paper we described the implementation of a Nitrogen III 
model ion into the NLTE atmosphere/spectrum synthesis code 
fastwind, to allow for future analyses of Nitrogen abundances in 
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O-starsO In particular, we concentrated on a re-investigation of 
the N in /L14634-4640-4642 emission line formation. Previous, 
seminal work by MH has suggested that the formation mecha- 
nism of these emission lines in 0((f)) and O(f) stars is domi- 
nated by two processes: the overpopulation of the 3d level via 
dielectronic recombination, and the strong drain of the 3p level 
(towards 2p 2 ) by means of the 'two electron transitions'. 

To account for the dielectronic recombination process in our 
new N in model, we tested two different implementations. Within 
the implicit method, the contribution of dielectronic recombina- 
tion is included in the photoionization cross-sections (provided 
from OPACITY project data), whereas the explicit method uses 
resonance-free photoionization cross-sections and directly ac- 
counts for the different stabilizing transitions from autoionizing 
levels. Both methods have been implemented into fastwind, and 
the equivalence of both methods has been shown. 

First tests were able to reproduce the results by MH quite 
well, after adapting our new atomic data to values as used by 
them. This result was achieved by applying similar assumptions, 
i.e., wind-free atmospheres and the same approximate treatment 
of line-blanketing. The triplet emission becomes even increased 
when using our new atomic data, mostly because of higher os- 
cillator strengths for the 'two electron transitions'. 

After switching to fully line-blanketed models (still wind- 
free), however, the situation changes dramatically. Then, the 
strong emission present in the 'pseudo line-blanketed' models is 
lost (consistent with results from the tlusty OSTAR2002 grid). 
We have investigated this mechanism in detail, and shown that 
it is rooted in the overall lower ionizing EUV-fluxes (not only 
until the Nm edge as assumed by MH), either directly or in- 
directly, and the strong collisional coupling of the ground state 
and the 2p 2 levels. In any case, these levels become much less 
underpopulated in fully line-blanketed models, thus preventing 
an effective drain of the 3p level. We also questioned the role 
of dielectronic recombination regarding the overpopulation of 
the 3d level. After setting the DR contribution to the 3d level to 
zero, almost no reaction on the triplet lines was found throughout 
the whole model grid. Both results lead to the conclusion that, 
under realistic conditions of line-blocking/-blanketing with so- 
lar background abundances, both the drain and the dielectronic 
recombination lose their key role as assigned to them by MH. 

This key role is now played by the stellar wind. Already MH 
suggested that the emission lines of Of sta rs (contrasted to 0((f)) 
and O(f) stars) might be formed due to fhe lSwingi d 19481) mech- 
anism, an overpopulation of the 3d level by pumping through the 
corresponding resonance line in an "extended atmosphere". This 
suggestion could not be proven though, since MH had not the 
tools to model such atmospheres. Nowadays, this is no longer a 
problem, and when we include the wind in our calculations (with 
mass-loss rates following the 'unclumped' Galactic WLR), it 
turns out that we obtain almost the identical emission as present 
in the original MH simulations (performed without wind and 
without realistic line-blanketing). By inspection of the net rates 
into and from the 3d level, we noticed how the wind induces the 
overpopulation of the 3d level via pumping from the ground state 
rather than by dielectronic recombination. A prerequisite of this 
process is that the wind-strength is large enough to display a sig- 
nificantly accelerating velocity field already in the photospheric 
formation region of the triplet lines, to allow the resonance lines 
to leave detailed balance. 



27 Corresponding Nn/Niv/Nv model ions will be presented in 
Paper II. 



The most important implication of our study is that under 
Galactic conditions DR plays only a secondary role, both for O- 
stars with "compact" and with "extended" atmospheres, whereas 
the key process is the overpopulation due to the resonance line in 
the presence of an accelerating velocity field. Note that without 
such a velocity field the resonance line loses its impact though. 
Consequently, it is to be expected that hydrostatic NLTE codes, 
such as tlusty and detail/surface, will not be suited to quanti- 
tatively synthesize the N in triplet lines (and, to a lesser extent, 
also N ni/i4097, because of their interrelation), unless the wind- 
strength is significantly below the Galactic WLR or the back- 
ground metallicit y, Z, is much lower than the Galactic one (see 
Heap et al.ll2006l who performed Nitrogen diagnostics for SMC 
O-stars by means of tlusty). Note also that our results have been 
derived for O-star conditions only, and should be valid as long as 
the N in emission lines, to a major extent, are formed in the pho- 
tosphere or in the transition region. For objects with significantly 
denser winds, e.g., WN-stars, our analysis would certainly need 
to be extended, since additional effects might be present or might 
even dominate. 

In order to check our new model atom and the predictive 
power of fastwind, we performed first comparisons with results 
from the alternative model atmosphere code cmfgen, for a small 
grid of O-type dwarfs and supergiants. For this objective we used 
a set of important N n and N in lines in the blue part of the visual 
spectrum. The overall agreement between both codes is mostly 
satisfactory, though some systematic deviations demand a fur- 
ther clarification in terms of a comparison with observations. 

Within the range 30,000 K < r eff < 35,000 K, however, 
some major discrepancies have been found. Here, cmfgen trig- 
gers the emission at - 4640 - 4642 earlier, i.e., at cooler 
temperatures, than calculations from fastwind. This effect could 
be traced down to the overlap of the N m resonance line (ac- 
tually, one of its fine-structure components, at 374.43 A) with a 
resonance line from O m, which is of similar strength throughout 
the O-star domain. We studied this effect by means of the cor- 
responding cmfgen calculations (since fastwind lacks a detailed 
Om model), and also by analytic considerations. 

As long as both resonance lines are optically thick, source 
function equality is achieved. Under most conditions, the 'cou- 
pled' source function lies in between the individual, 'decoupled' 
ones from N m and O m. In the critical temperature range now, 
the decoupled O ra source function is predicted to be much larger 
than the N m one (due to substantial cascades from upper levels), 
and leads, after being coupled with Nm, to significant values 
for the combined source function, beyond the decoupled Nm 
case. Thus, the important 3d level becomes more pumped, and 
the triplet emission occurs at cooler temperatures than without 
coupling. When the coupling is neglected in cmfgen, the pre- 
dicted triplet emission vanishes and the profiles become similar 
to those from fastwind. 

It might be suspected that the impact of the O m resonance 
line overlap introduces an additional parameter to be known 
when the Nitrogen abundance shall be determined via the triplet 
lines, namely the Oxygen abundance. As we have shown by our 
analytic considerations, however, this parameter remains rather 
unimportant as long as both resonance lines are optically thick, 
which is true in the interesting parameter range. We have tested 
this prediction by lowering the O m oscillator strength in cmfgen 
by a large factor (50), and found no difference in the coupled 
source function and triplet emission. 



except for certain differences in the EUV-fluxes etc. 
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Summarizing, not only the N m resonance line itself is re- 
sponsible for the (strong) triplet emission, but also the overlap- 
ping O in resonance line, at least at later O-types. Indeed, for one 
corresponding object studied in Paper II, Nl 1-029 (09.71b), we 
encountered the problem that the observed N in triplet could not 
be reproduced by fastwind, though with cmfgen, and we inter- 
preted this problem as due to the resonance line overlap. Insofar, 
we need to incorporate this process into fastwind if we aim at 
deriving Nitrogen abundances at cooler temperatures from the 
triplet lines alone. Fortunately, the contamination of the other 
optical lines by this process remains weak, so that these lines 
can be used already now, and the Nitrogen abundance determi- 
nation by fastwind is not hindered. 

Our study implies two important consequences that need to be 
tested in future comparisons with observations. 

(i) Since the efficiency of DR and 'two electron' drain 
is strongly dependent on the degree of line-blanketing/- 
blocking, we predict that in a metal-poor environment (e.g., 
in the Small Magellanic Cloud with Z/Z « 0.2) the emis- 
sion becomes stronger again, due to less EUV line-blocking. 
On the other hand, in such a low-Z environment also the 
wind-strengths and the base-line Nitrogen abundance be- 
come lower, and the combined effects need to be investigated 
in detail. 

(ii) As outlined above, the triplet lines from O-type stars are of 
photospheric origin and depend, via the N m (and O m) res- 
onance lines, on the actual wind-strength (independent of 
clumping and X-ray properties), determining the onset of 
the accelerating velocity field. Thus, their emission strengths 
might be used to constrain the stellar mass-loss rate (af- 
ter the line-formation has been shown to work reliably), 
if the Nitrogen abundance can be derived independently 
from other lines. In particular, 'weak-w i nded' stars (e.g. , 
Bouret etafl 120031: iMartins et all 120051* IPuTs et all l2008t 



Marcolino et al. 



120091: iNaiarro et all 1201 lh have a rather 
weak wind for their luminosity, being a factor of 10 to 
100 thinner than predicted/observed for their counterparts 
with 'normal' winds. Thus, it might be suspected that weak- 
winded stars display much less emission than stars with 'nor- 
mal' winds of similar type. This requires, that a) there are 
weak-winded stars also at intermediate O-types, and b) that 
the winds of their 'normal' counterparts are not strongly 
clumped, which would diminish also their emission because 
of lower-than-thought mass-loss rates. 
Using the N m triplet as an independent mass-loss diagnos- 
tics would be somewhat simi lar to the correspo nding ap- 
plication of_the NIR Br„-line dNajarro et alJl20lH see also 
iPuls etal] [20081. but with the advantage that the emission 
strength of the triplet is a rather monotonic function of M, 
whereas Br tJ , changes its behaviour from weak to normal 
winds considerably. 
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Appendix A: Dielectronic Recombination: 
Implementation to fastwind 

In the present work we implemented dielectronic recombination 
into fastwind, which, so far, could not (or only approximately) 
deal with this process. To this end, new rates (both for the di- 
electronic recombination and for the inverse process) had to be 
inserted into the system of the rate equations. 



A. 1. Explicit method 

To calculate these rates for the 'explicit method' (see Sect. I3.2I>. 
we fol low the formulation as provided bv lNussbaumer & StorevI 
d!983l) . In compact notation, the dielectronic recombination for 
an element X and charge m + 1 proceeds via 



X" 



XI 



X™+hv 



(A.l) 



where p is a parent state of the m + 1 times ionized element X, 
a is an autoionizing state and b is a bound state. We denote the 
initial state of expression lA.il composed of the recombining ion 
and the free electron, as a continuum state c. We refer to the first 
process as dielectronic capture and to its inverse as autoioniza- 
tion. In general, dielectronic captures and autoionizations link 

state a to a large numb er of continuum states c . 

As a final resultj^ Nussbaumer & Storey (1983) could ex- 
press the dielectronic recombination coefficient between au- 
toionizing state a and bound state j, aoR(a j), as 



n e N™ +l a DR (aj) = b a {n™yA*, 



(A.2) 



with A*, the corresponding radiative transition probability for 
the stabilizing transition (Einstein coefficient for spontaneous 
emission, corrections for induced emission will be applied be- 
low), total ion density A^" +I (element k), («™)* the LTE popu- 
lation of state a with respect to the actual (NLTE) ground-state 
population of the next higher ion and the actual electron density, 
and b a the related NLTE departure coefficient. 

The last quantity can be expressed in terms of (i) the au- 
toionization coefficient(s), A a ac , between state a and all possible 
compound-states c that can form a, (ii) the radiative transition 
probabilities, A\ between state a and all possible bound and au- 
toionizing states with lower energy i to which state a can decay, 
and (iii) the departure coefficients of the contributing parent lev- 
els, bp (here with respect to the ground-state of the same ion, 



(A.3) 



Usually, the autoionizing probabilities for state a are much larger 
than the radiative probabilities for decay, and often there is only 
one parent level, namely the ground-state of ion m+l, n p = n™ +l , 
i.e., bp = 1. Under these conditions (which are similar for ex- 
cited parent levels assumed to be in LTE with respect to the 
ground level), b a —> 1, and the dielectronic rate depends only on 
the LTE population of state a and the radiative transition prob- 
ability A R .. All dependencies on the autoionization probabilities 
have 'vanished', and we need only the value of («'")* that can be 
derived from the Saha-equation and the ground-state population 
of ion m + 1 , without including the autoionizing levels into the 
rate equations. 



See also Mihalas 1978 for a simplified derivation. 



In stellar atmospheres, one needs (in addition to the sponta- 
neous emission to level j) to account for stimulated emission as 
well, i.e., 



Kj + BajJ 



1 + 



2hv 3 



(A.4) 



with B a j the Einstein coefficient for stimulated emission and J 
the scattering integral (profile weighted, frequency integrated 
mean intensity) for the stabilizing transition. Since the important 
resonances are broad, the scattering integrals might be replaced 
by the mean intensities, J v , of the pseudo -continuum (i.e., in- 
cluding all background opacities/emissivities) at the frequency 
of the stabilizing line. 

Finally, we can define the total dielectronic rate to level j 
from any possible autoionizing state a', 



n e N™ +1 X ataztfj) * ne< +1 -CJ-^x (A.5) 

*f 8k 

x J^op(-^/ir)A; j (l + y-rJil). 

i U 

which might need to be augmented by departure coefficients b a , 
inside the rhs sum if the parent levels are not the ground-state or 
not in LTE with respect to the ground-state. The inverse upward 
rate involves only the line process(es), 



< y B, a iJij = n>- y -^-A* .^/ y> 



Si 



(A.6) 



with Bjj the Einstein-coefficient(s) for absorption. It is easy to 
show that for LTE conditions (n'J - (n™)* and Planckian radi- 
ation in the lines) the upward and downward rates cancel each 
other exactly, as required for thermalization. 

Though we have fol lowed here the derivation by 
Nussbaumer & StorevI (119831) . our final results for the downward 
and upward rates are ident ical with the formulation as used by 
Mihal as & Hummer! ( 1 1973b in their investigation. 

These rates have been implemented into fastwind and are 
used whenever the 'explicit method' is applied. The only input 
parameters that need to be specified in the atomic data input file 
are the transition frequencies and the oscillator-strengths for the 
stabilizing lines, fj a , which relate to the products g fl A R . via 



%n l e 



,2„2 



aj m e c 3 aj 



(A.7) 



For convenience and for consistency with the formulation of 
'normal' recombination rates (see Sect. IA.2l >. the quantity («"')* 
is expressed in terms of the LTE -population (with respect to ac- 
tual ionization conditions) of the lower, stabilizing level (n" 1 )* 



«)* = (njy^exp(-hv aj /kT), 
8j 



(A.8) 



so that the downward rate (for a specific autoionizing level a) 
can be expressed as 



2„2 



J m e c J 1 



x ( ! + 2^4 



(A.9) 



Summation over all contributing autoionizing states yields the 
final rate. 
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In this explicit method, the rates for dielectronic recombi- 
nations and inverse processes are calculated in a separate step 
and then added to the rates involving resonance-free photoion- 
ization cross-sections alone. In our model ion (Sect. 0), we use 
such cross-sections defined in terms of the Seaton approximation 
(Eq. |H, which, together with the data for the stabilizing transi- 
tions, have been taken from the WM-basic atomic database. Note 
that we consider processes both from/to ground-states as well as 
from/to excited states within ion m + 1, so far on the assumption 
that the autoionizing levels are in LTE (i.e., without including 
these levels into the model atom and rate equations). 

A.2. Implicit method 

Within the implicit method, the DR contribution is al- 
ready contained within the 'conventional' recombination rates, 
(nj/rikfRkj, to yield a total number of recombinations 



n k 



(A. 10) 



As usual, n/t is the actual population of the recombining ion in 
state k and the LTE population ratio of the considered 

level to which the process recombines (either directly or via the 
intermediate doubly excited state). Rkj is defined as 



R kj = 4tt 



2hv i 



,-hv/kT 



dv 



(A.ll) 



with mean intensity J v and total photoionization cross-section 
(including resonances), acjk(y). 

In the following, we show that this formulation is consistent 
with the rates derived in Sect. I A. 1 1 We split the cross-section 
into a resonance-free contribution, and a contribution from the 
resonances, 



(A. 12) 



The total recombination rate is then the sum of direct and dielec- 
tronic recombination, 



nfRkj 
with 



< S (v) / 2hv 



r 

Jres hv 



-hy/kT dy _ 



(A.13) 



(A. 14) 



The resonances are narrow enough so that most of the frequency 
dependent quantities can be drawn in front of the integral, eval- 
uated at the average position of the resonances i, 



4 n 2hvl _ 
1 hvu c- 



hVij/kT 



(A.15) 



The integral over the cross-sections of the resonances corre- 
sponds to the cross-section of the stabilizing transitions, 



An f 

7 ^ 

nv U Jres(i) 



4n ne 
dv = Bp = fy, 

nvij m e c 



(A. 16) 



with Einstein-coefficient Bj, and oscillator-strength fy. Like- 
wise, 



«7res(i) 



aJJ v dv oc fyJij. 



Then, indeed, we recover the result from Eq. IA.9I (explicit for- 
mulation), 



v-i 87T e~ 9 



2, ,2 



1 + 



2hv 



3 J U 



(A. 17) 



if, as outlined in Sect. IA. II the autoionizing levels are in or close 
to LTE. Note that this is a necessary condition for the implicit 
method to yield reliable results^ otherwise one has to use ex- 
clusively the explicit approach and to include the autoionizing 
levels into the model atom and rate equations. 

The proof that the rates for the inverse photoionization pro- 
cess, calculated either by the implicit method (using total pho- 
toionization cross-sections) or via rates from resonance-free 
cross-sections plus rates involving the excitation of the second 
electron, are consistent, is analogous and omitted here. 

A.3. Consistency check 

To check the consistency of our implementation of implicit and 
explicit DR treatment, we carried out the following test. First, 
to ensure consistent data, we derived the oscillator strength cor- 
responding to the wide PEC resonance in the total photoioniza- 
tion cross-sectio n (Fig. E ll by integrating over the cross-section 
and applying Eq. IA.16PI The resulting value (/ = 0.45, which 
is somewhat smaller than the data provided by the WM-basic 
database, / = 0.60) was then used within the explicit method, 
at the original wavelength (which coincides with the position 
of the resonance). As displayed in Fig. I A. 1 1 for the case of 
model 'T3740' with 'pseudo line-blanketing' (see Sect. 15. IK 
both methods result in very similar line profiles for the N m emis- 
sion triplet, proving their consistency. Figure I A. 1 1 also displays 
the strong reaction of the N m triplet when the oscillator strength 
is either increased or diminished by a factor of two. 



30 which has been used to calculate the total cross-sections as provided 
by, e.g., the OPACITY Project. 

31 Of course, the underlying contribution responsible for direct ion- 
ization needs to be subtracted. 

NIII4634-4640-4642 




4632 4634 



4636 4638 4640 
wavelength (A) 



4642 4644 



Fig.A.l. Consistency of implicit and explicit method. 
Comparison of Nm /U4634 - 4640 - 4642 profiles from 
model 'T3740', for calculations using a different treatment 
of DR to the 3d level. Implicit method (solid/black), explicit 
method (dotted/red), explicit method with oscillator strength for 
the stabilizing transition increased (dashed/blue) and decreased 
(dashed-dotted/black) by a factor of two. 
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Fig. B.l. Grotrian diagrams for the N in doublet (left) and quartet (right) system. Level numbers refer to Table IB. fl Levels #8, 10, 
1 1 refer to the 3s, 3p and 3d levels involved in the formation of N m AA4634 - 4640 - 4642, and # 3, 4, 5 refer to levels 2p 2 ( 2 D, 2 S, 
2 P) that drain level 3p via 'two electron transitions' . Important optical transitions are indicated by green lines and numbers referring 
to Table |4]in the main section. 



Appendix B: Details on the Nm model ion 

Table lB. fl displavs the electronic configurations and term desig- 
nations of our Nm ionic model, and Fig. IB. li the corresponding 
Grotrian diagrams for the doublet and the quartet systems, indi- 
cating important diagnostic transitions in the optical (cf. Sect.|6]i. 



all models from the grid presented in Table [3] For a discussion, 
see Sect. [6] 



Table B.l. Electronic configurations and term designations of 
our N in model ion. The level numbers correspond to the entries 
in the Grotrian diagrams, Fig. IB. II 



# 


Configuration 


Desig. 


# 


Configuration 


Desig. 


1 


Is 2 


2s 2 ('S)2p 


2p i P ,) 
2p 2 4 P 


22 


Is 2 


2s 2p( J P u ) 3p 


3p' 2 D 


2 


Is 2 


2s 2p 2 


23 


Is 2 


2s 2p( 3 P°) 3p 


3p' 2 S 


3 


Is 2 


2s 2p 2 
2s 2p 2 


2p 22 D 


24 


Is 2 


2s 2p( 3 P°) 3d 


3d' 4 F° 


4 


Is 2 


2p 22 S 


25 


Is 2 


2s 2p( 3 P°) 3d 


3d' 4 D° 


5 


Is 2 


2s 2p 2 


2p 2 2 P 


26 


Is 2 


2s 2 ('S) 5s 


5s 2 S 


6 


Is 2 


2 P 3 


2p 3 4 S° 


27 


Is 2 


2s 2p( 3 P°) 3d 


3d' 2 D° 


7 


Is 2 


2 P 3 


2p 3 2 D° 


28 


Is 2 


2s 2p( 3 P°) 3d 


3d' 4 P° 


8 


Is 2 


2s 2 ('S) 3s 


3s 2 S 


29 


Is 2 


2s 2 ('S) 5p 


5p 2 P° 


9 


Is 2 


2 P 3 


2p 3 2 P° 


30 


Is 2 


2s 2p( 3 P°) 3d 
2s 2 ('S) 5d 


3d' 2 F° 


10 


Is 2 


2s 2 ('S) 3p 


3p 2 P° 


31 


Is 2 


5d 2 D 


11 


Is 2 


2s 2 ('S) 3d 


3d 2 D 


32 


Is 2 


2s 2p( 3 P°) 3d 


3d' 2 P° 


12 


Is 2 


2s 2p( 3 P°) 3s 


3s' 4 P° 


33 


Is 2 


2s 2 ('S) 5f 


5f 2 F° 


13 


Is 2 


2s 2p( 3 P°) 3s 


3s' 2 P° 


34 


Is 2 


2s 2 ('S) 5g 


5g 2 G 


14 


Is 2 


2s 2 ('S) 4s 


4s 2 S 


35 


Is 2 


2s 2 ('S)6p 


6p 2 P° 


15 


Is 2 


2s 2p( 3 P°) 3p 


3p' 2 P 


36 


Is 2 


2s 2 ('S) 6s 


6s 2 S 


16 


Is 2 


2s 2p( 3 P°) 3p 


3p' 4 D 


37 


Is 2 


2s 2p('P°) 3s 


3s" 2 P° 


17 


Is 2 


2s 2 ('S) 4p 


4p 2p0 


38 


Is 2 


2s 2 ('S) 6d 


6d 2 D 


18 


Is 2 


2s 2p( 3 P°) 3p 


3p' 4 S 


39 


Is 2 


2s 2 ('S) 6f 


6f 2 F° 


19 


Is 2 


2s 2p( 3 P°) 3p 


3p' 4 P 


40 


Is 2 


2s 2 ('S) 6g 


6g 2 G 


20 


Is 2 


2s 2 ('S) 4d 


4d 2 D 


41 


Is 2 


2s 2p( 3 P°) 4s 


4s' 4 P° 


21 


Is 2 


2s 2 ('S)4f 


4f 2 F° 











Appendix C: Comparison of line profiles with 
results from cmfgen 

In this appendix, we display a comparison of N in (and partly 
N n) line profiles from fastwind (black) and cmfgen (green), for 
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Fig. C.l. Comparison of Nn line profiles from present work 
(black) and cmfgen (green), for model dlQv (for parameters, see 
Table©. 
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Fig. C.l. As Fig. [CTTl but for model d8v. 
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Fig. C.3. As Fig.lCHl but for model slOa. 



Fig. C.4. As Fig. |CT] but for model s8a. 
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Fig. C.5. Comparison of Nm line profiles from present work 
(black) and cmfgen (green), for model dlOv. 
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Fig. C.6. As Fig. |C5] but for model d8v. 
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Fig. C.7. As Fig.|C3] but for model d6v. 



Fig. C.10. As Fig. |C3] but for model slOa. 
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Fig. C.8. As Fig.[C5] but for model d4v. 
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Fig. C.ll. As Fig. |C3] but for model s8a. 
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Fig. C.9. As Fig. |C31 but for model d2v. 
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Fig. C.12. As Fig.lC31 but for model s6a. 
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Fig. C.13. As Fig. |C3J but for model s4a. 
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Fig. C.14. As Fig.|C5] but for model s2a. 



